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^ I Abstract 
P 

The production of dileptons is studied within a hadronic transport model. 
We investigate the sensitivity of the dilepton spectra to the initial configuration 
of the hadronic phase in a ultrarelativistic heavy ion collision. Possible in 
medium correction due to the modifications of pions and the pion form factor 
in a hadronic g discussed. 
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1 Introduction 



One of the major goals of the ultrarelativstic heavy ion program is to study the restoration 
of chiral symmetry at high temperatures and densities and possibly to create and identify 
a new form of matter the so called Quark Gluon Plasma. To this end the measurement 
of electromagnetic probes such as photons and dileptons have continually attracted great 
interest. Contrary to hadronic observables these weakly interacting probes are not af- 
fected by final state interactions and thus may provide information about the early, hot 
stage of the reaction. While at very high invariant masses, above that of the J/ty the 
dilepton yield is essentially dominated by the Drell-Yan production from the hard nucleon 
nucleon collisions, it has been suggested that in the mass region between the $ and the 
J/^f, dileptons originating from the QGP can be observed (for an overview see e.g. ref. 
0). Of particular interest in the high invariant mass region has been the study of the 
suppression of the J/\P, because of its possible implications about deconfinement. At low 
invariant masses, below 1 GeV the dilepton spectrum is expected to be dominated by the 
hadronic phase, where processes such as hadronic decays and pion annihilation contribute. 
This region is of particular interest, because possible in medium modifications of hadron 
masses such as that of vector mesons may be observable in the dilepton invariant mass 
spectrum 0, [3], [7j . This mass region is also accessible at lower bombarding energies (BE- 
VALAC, SIS) where modifications of hadron masses due to finite nuclear density rather 
than temperature can be studied [[8]. [| 

Recently interest |]TT , |l~2] , p2| , [13], 1|| in the low mass region has been sparked by the 



experimental data from the CERES collaboration ||, which show a considerable enhance- 
ment at invariant masses of about 400 MeV for SPS-energy S+Au heavy ion collisions as 
compared to proton-proton as well as proton-nucleus collisions. A similar enhancement 
has also been seen by the HELIOS collaboration at more forward rapidities [jl0|1 . It has 



been suggested |LJ], [12] that this enhancement may be due to the dropping of the vector 
meson masses at finite temperature and densities, as proposed in [[15], [L6|. However, also 
more conservative effects, such as a modified in medium pion dispersion relation may give 
sufficient enhancement ||, This later effect has first been suggested in the context of 
BEVALAC energy heavy ion collisions ]T7|, [IS, |T^] where the pion dispersion relation is 



modified due to the interactions with nucleons forming deltas. At SPS energies, one would 
expect that the major modification should come from interactions among pions, mostly 
in the the isovector - p-wave channel, which is dominated by the s-channel rho-meson 
resonance [pC, 2 1 1 . The effect of the nucleons, however, may not be negligible [22]. 

However, before definite conclusions about possible in medium effects of any kind can 
be drawn, a conventional explanation of the observed enhancement has to be ruled out. 
The cocktail used by the CERES collaboration to estimate the conventional background 
includes only the decay of hadrons with a relative abundance taken from proton proton 
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data and scaled with the number of charged particles. In a heavy ion collision several 
things could be different from these assumptions due to multiple elastic and inelastic 
scatterings. First of all the relative abundance could be altered due to creation an de- 
struction of higher lying meson resonances such the p or the a%. Secondly, there could be 
multiple 7r-7r annihilation processes leading to additional dilepton yields. Furthermore, 
initially , after formation, the hadrons not necessarily have to be in thermal equilibrium. 
One rather could imagine that a great deal of the thermalization observed in the final 
(freeze out) state |23| actually takes place within the hadronic phase. This equilibration, 
which involves elastic and inelastic collisions, then may result in a quite different dilepton 
yield as, for instance, obtained from an equilibrated hadron gas. In other words, a careful 
and systematic investigation of the dilepton production may reveal information about the 
equilibration properties of the hadronic phase. 

It appears rather unlikely, on the other hand, that the equilibration process can be 
inferred from hadronic observables only. Provided that a certain degree of equilibration 
has been obtained in the final state, as suggested by a recent analysis of AGS as well 
as SPS data ||24 , |23| , it is of course impossible to extract in which way the system has 
reached this point. Penetrating probes, however, are not only sensitive to the final state 
but represent an integral over the entire reaction history and thus are much better suited 
to constrain our understanding of equilibration dynamics. In other words, the fact the 
event generators such as RQMD p5) or VENUS p6f are able to reproduce most hadronic 
observables not necessarily implies that they simulate the reaction dynamics properly, 
since there are infinitely many trajectories towards thermal equilibrium. 

The purpose of this article is to investigate the sensitivity of the dilepton yield to 
different initial conditions in the hadronic phase and to see to which extent the CERES and 
HELIOS data can be accounted for by non equilibrium processes in the hadronic phase. 
Our strategy therefore is to parameterize the initial phase space and the abundances 
of the hadrons subject to the constraint, that the final hadronic spectra and rapidity 
distributions are reproduced. This will lead to a set of initial conditions for which we 
then study the contribution to the dilepton data. 

This paper is organized as follows: In the first section we will describe the hadronic 
transport model. In the second section we will describe the initialization of the hadronic 
phase space. In the third section the production of dileptons within out transport model 
will explained. We then will turn to the discussion of possible medium effect before we 
present the results and our prediction for Pb+Pb collisions. 
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2 Transport Model 



In this section we will describe the hadronic transport model used to study the dilepton 
production. In our calculation we include pions, p- u>- and a\ mesons as well as nucleons 
and deltas for the baryons. We ignore higher baryon resonances at this time but we 
will present an estimate of their contribution to the dilepton spectrum below. Since the 
fraction of baryons in the relevant rapidity region is small compared to that of the pions, 
we do not expect that they contribute significantly to the equilibration process. We have 
also left out the strange mesons for the same reason. The contribution of the eta meson 



to the dilepton yield is calculated at the end using the measured eta/pi ratio [27j. 

The interaction between nucleons, deltas and pions is done in the standard fashion as 
reported e.g. in p8[ ^] and we refer the reader to these references for details. 



All hadronic resonances are propagated explicitly. For instance the resonant isovector 
p-wave 7T — 7T scattering in the p channel is split into the formation, the propagation and 
the decay of a p-meson. As a consequence the mass of these resonances are not restricted 
to their central value but may assume any value as determined by the production process 
of by the initialization procedure (see below). To be specific, for the formation of the p 
from 7i — 7i collisions, we use an isospin averaged Breit-Wigner cross section 

a 7T7r ^ p {M) 



p 2 (m 2 - m 2 ) 2 + M 2 r 2 (M) v ' 

Here p is the center of mass momentum and M the invariant mass of the pair. Using the 



standard Ti-p coupling [311 the width of the p is given by 



r P _(M)=r 0iP ^ (2) 

where p is the decay momentum. T 0>p = 150 MeV and p are the width and the decay 
momentum at peak value (M = m p ). 

The formation and decay of the a x is a little more complicated since it involves a 
hadronic resonance, the p, in the initial or final state, respectively. Thus, the a\ does 
not necessarily decay into a p of mass m p = 770 MeV but also into states with smaller 
and larger masses. As a starting point, we use the formula of reference [32] for the decay 
width of the a x into pion and p 

r« (M, m p ) = Q(M 2 - ml - m 2 ) 2 + mj(m 2 + p 2 )) (3) 



«i 



Notice, that in the above formula the width depends also on the mass of the outgoing p. 
Therefore, if we want to allow for the decay into p-mesons with masses different from the 
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central value the total ai-width should be given by the integral of the above formula over 
all kinematically allowed p-masses weighted with the p mass distribution f p (m). Hence 

roc 

ltf(M)=/ f p (m p )T ai (M,m p ) (4) 

with 

2 m 2 T p (m) 
tp[m) = n (m 2 - m 2 ) 2 + m 2 F 2 (5) 

We thus have to adjust the coupling G in eq. (||) such that the total width Y l °^ agrees 
with the experimental value of 400 MeV. We find G aip7T = 18.1 GeV -2 . The partial width 
for the decay into a p of a specific mass is then given by 

dT 



dm p 



f p (m p )T ai (M,m p ) (6) 



This then gives the distribution from which the mass of the outgoing p has to be picked. 
Finally, the isospin averaged cross section for the formation of an a\ in a n — p collision 
is given by 

3vr M*T ai (M,m p )T%(M) 
cr^aA^) 4p2 (M2 _ m2 j 2 + M 2 V m^ M) V) 

where T ai (M, m p ) is given by eq. ([|) and r^°*(M) is given by eq. fl5|). The advantage of 
this somewhat complicated prescription is that it allows to include the finite width of the 
p and, at the same time, preserves detailed balance. 

For the u meson we adopt a similar strategy. We model the dominant decay into 
three pions as a two step process via an intermediate 7r p state, as suggested by chiral 



models including vector mesons [|33|] . Thus in our approach the u decays into a pion and 
a p-meson the mass of which is below the central value. As explained previously in our 
transport model these p's are treated as real particles. To preserve detailed balance we 
also allow for the formation of the uj by it p-fusion. 



Following ref [33] we use a p-wave form for the decay width of the u 



r w (m p ) = Cq 3 (8) 

where q is the decay-momentum of the omega in its rest frame. The constant of propor- 
tionality, C, is then determined such that the total width, integrated of all possible p - 
masses with the proper weighting, corresponds to the observed experimental value for the 
decay into three pions 



^tot 



C f p (m p )q 3 dm p = 7AMeV (9) 
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which gives a value of C = 48 GeV -2 . 

Following the same arguments as in case of the a% the isospin averaged cross section 
for the formation of the u is then given by 

g (M, m ,) = il r " (m ;> r " , (10) 
v ' p> 9q 2 (M 2 - m 2 ) + M 2 r*f v ' 

Naturally, due to the small width of the u> very few new u are being formed in the 
course of the expanding hadronic system. Therefore, the resulting dilepton spectrum is 
only slightly changed as compared to calculations which do not take the formation of the 
omega into account JTlJ . 



Finally, for the non resonant scatterings we assume a constant cross section of 20 mb 
except for that among nucleons and Deltas, where a parameterization of the measured 
nucleon-nucleon cross section is being used. 



3 Initialization of initial Phase space 

As already outlined in the introduction, in our transport model we will not describe the 
entire collision history, starting from the initial colliding nuclei. We rather want to ask 
the question about the sensitivity of the dilepton production to the hadronic phase-space 
configuration created in such a collision. Thus the initial stage in our description is that 
of a hadronic gas with a possible non equilibrium configuration in momentum space. We 
also do not impose chemical equilibration. We, however, require that measured hadronic 
spectra and rapidity distributions are reproduced at the end of the evolution. As we shall 
see this leaves quite a number of possibilities to initialize the system. 

3.1 Configuration Space: 

We will assume that initially all particles are distributed uniformly within a cylinder of 
radius Rq and longitudinal extent 2Zj. This somewhat simplified prescription may be 
improved by e.g. using a Gaussian distribution in the longitudinal direction, which we 
do not consider here. However, we should point out that the introduction of a formation 
time, discussed below, will effectively introduce some smearing of the particle density in 
the longitudinal direction. Since we plan to compare with the CERES data for S + Au 
collisions, we will assume an initial transverse radius of Rq = 3.5 fm. The longitudinal 
extent will be varied and the specific choices will be given in the results section (||). 
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3.2 Momentum Space: 



Momenta in the transverse direction are distributed according to a two dimensional Bose- 
or Fermi- distribution with the possibility of non vanishing chemical potential. The chem- 
ical potential of the p is then twice that of the pion and those of the u and a% three times 
the pionic one. 

For the longitudinal momenta we assume that the rapidities of the particles are dis- 
tributed according to a Gaussian distribution, 



which allows to take into account possible correlations between rapidity and longitudinal 
position z. To some extent, correlations between the longitudinal position and the rapidity 
are already generated by the introduction of a formation time as we will discuss below. 

For the hadronic resonances the initial distribution of the masses has to be specified as 
well. In principle there are two possibilities. First, the masses are distributed according 
to the respective Breit Wigner distribution fhadron( m )- Second, assuming that thermal 
equilibration has been achieved prior to hadronization, the masses should be distributed 
according to thermally weighted distribution, which slightly favors lower mass states. 
Assuming boost invariance on gets 



where T is the transverse temperature, K 2 is the modified Bessel-function, and N is a 
normalization factor. Notice, that one obtains the same mass distribution in case of a 
Bjorken type fire cylinder. In this work we will only consider the second possibility ([12]). 

3.3 Formation Time: 

A realistic description of the initial hadronic fireball should include a formation time. 
Particles with higher rapidities will be created at later times in a given frame of reference. 
In our approach we assume that particles will propagate without interaction for times 
smaller than their formation time 



where y is the rapidity of the particle under consideration and r denotes the proper 
formation time. This formation time clearly introduces a smearing of the longitudinal 
distribution in configuration space in the sense that, at the time particles will be allowed to 
interact, the fastest ones have already moved a considerable distance 5z ~ tf orrnation c. For 
the same reason the formation time leads to a rapidity longitudinal-position correlation. 






(12) 



t formation TCOsll(?/) 



(13) 
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3.4 Particle Abundances: 



In order to reproduce the measured rapidity spectrum of protons by the NA35 collabora- 
tion 0] we assume that our hadronic fireball contains 60 baryons. The relative weighting 



is then determined by the temperature assuming thermal and chemical equilibrium. At 
the temperatures considered here we will have about as many nucleon as deltas. 

For the mesons we will either use the relative abundances as determined from proton 
proton collisions, i.e Np/N^ = 0.2, N^/Np = 1/3. Since the number of at is not given by 
experiment we assume N ai /N u = 0.5. 

Alternatively, we will assume initial chemical equilibrium and will determine the rel- 
ative abundances assuming a boost invariant momentum distribution. In this case the 
abundances of course depend on the initial temperature and will therefore be given in the 
appropriate context in the results section @. 

3.5 Flow: 

Several analyses of available SPS data suggest that there is substantial radial flow, which 
leads to larger apparent temperatures [^, |23|]. We, therefore, allow for the possibility 



that the hadrons have some transverse flow in the initial state. For simplicity we assume 
a constant flow velocity independent of the radial distance of the particle. To generate the 
flow, we first distribute the particles thermally in the transverse direction, then we boost 
them according with the flow velocity before we finally boost them in the longitudinal 
direction as given by the initial rapidity distribution. This sequence of boost is identical to 
that employed by Schnedermann et al. []3~5|| . Allowing for transverse flow reduces the initial 
temperature needed to reproduce the final transverse momentum spectra. We should 
point out, however, that in our transport model a certain amount of flow is generated as 
result of particle collisions. Therefore, the initial flow velocities used here should not be 
compared directly with results from the analysis of hadronic spectra. 



4 Dilepton production 

Having specified the initial phase space distribution as well as transport model we finally 
turn to the production of dileptons. As far as dilepton production through pion annihi- 
lation is concerned there are essentially two ways to proceed in a transport model, which 
are in principal equivalent. One is to directly calculate the dilepton yield from the pion 
collisions and the other one is to obtain the dileptons from the decay of the p which has 
been formed in the collisions. Both are equivalent as has been shown in |7|, [IT] . However, 



we found that the second prescription is extremely sensitive to statistical fluctuations for 
invariant masses smaller that 400 MeV, as we will explain below. We, thus, resort to the 
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first prescription, namely the production of dilepton directly from the tt tt collisions. We 
use the well known cross section 



2 47T a r Am w m 4 p 
a n+ ^ e+e -{M) - gYJp V ~ M 2 (M 2 - mj) 2 + M 2 T p (M) 2 { ] 

where the factor | arises from isospin averaging. 

All the Dalitz decays are done explicitly. Therefore, in order to avoid double counting, 
only the p mesons present in the initial state are allowed to decay into dileptons. All 
secondary p mesons come from processes such as pion annihilation or a% decay, where 
we take the contribution to the dilepton yield into account explicitly, including vector 
dominance, so that those p mesons should not decay into dileptons. 

We furthermore include the direct decay of the uj into dileptons as well as the Dalitz- 
decay of the 77 meson. As already pointed out previously, we do not explicitly propagate 
the rj, but include it at the end of the calculation. To that end we assume that the 
etas have the same rapidity distribution as the pions. The relative abundance and the 
transverse momentum spectrum we obtain from the measured 77/710 ratios as a function 



of pt provided by the WA80 collaboration [27 



In order to improve statistics we treat the dilepton decay 'perturbatively', i.e. at 
each time step the particles are allow to 'radiate' dileptons, without being removed. The 
contribution to the final spectrum is then given by the partial width multiplied by the 
time step size T e + e - dt. At the end of the calculation all remaining hadrons are forced to 
decay and their contribution to the dilepton spectrum is given by the respective branching 
ratios. In case of the pion annihilation, any ix tt collision contributes to the dilepton yield 
with the ratio of dilepton production cross section over total cross section. 

For a comparison with the CERES data, cuts on the lepton momenta have to be 
performed. In order to improve our statistics further we sample the decay of each virtual 
photon several (ten) times. 

For the Dalitz decays of the 77 and ui we use the formulae of Landsberg [p6f including 



form factors. For the form factors we use the values given in [36 



f v {m) = (l-2.4m 2 ) 2 (15) 
Urn) = (l.-m 2 /0.72 2 ) 2 (16) 

where the invariant mass m is to be given in GeV. 

For the a\ we also use the Landsberg formula with a vector dominance form factor 
including the width of the p 

4 

777 

= (M 2 - 777,4)2 + M 2 r p (M) 2 ^ 
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We have compared this formula with an explicit calculation of the Dalitz decay of the a\ 
using the Lagrangian of p2] an found very good agreement. 



For the partial decay width of the p we use ]5T 



771 

T po _ >e+e -(M)=T OJ £ (18) 

where T = 6.77 keV is the measured width at peak value. Notice, that for small invariant 
masses the branching ratio for a p decaying into dileptons increases rapidly since the total 
width, which is dominated by the p-wave decay into pion becomes small. As a consequence 
if an initial p with small invariant mass happens to survive until the end of the calculation 
it will decay with a very large branching ratio. We find that for 400 runs the average 
number of p with small invariant masses surviving until the end (30fm/c) is less than one 
and thus the statistical fluctuations are tremendous. By changing the random number 
seed the results for invariant masses smaller that 400 MeV fluctuated by more than an 
order of magnitude. Since this mass region is identical to the one where the CERES 
collaboration reports a strong enhancement, these uncertainties are unacceptable. We 
should also point out that we have observed similar fluctuations when we calculated the 
pion annihilation in a two step process, i.e. via the formation and the decay of the p meson. 
To avoid those we resorted to the direct production of dilepton from pion collisions as 
described above. In order to circumvent the statistical fluctuations from the p mesons 
in the initial state we calculate the contribution to the dilepton spectrum from p mesons 
with masses smaller than 500 MeV directly, i.e. not by Monte Carlo methods. To this end 
we calculate the number N p (M) of p meson in a given mass interval M + AM from the 
mass and phase space distribution. The dilepton yield is then obtained by multiplying 
with the branching ratio 

dNe+e ~ - N p {M) Vp ^ e+e ~ (19) 



dM 



tot 



Since N P (M) is proportional to T to t the above distribution is smooth even if T to t becomes 
small. By calculating the decay of these small mass p mesons directly we ignore possible 
absorption effects and thus somewhat overestimate the yield at low invariant masses. 
Naturally, in order not to double count, all initial p mesons with masses below 500 MeV 
are not allowed to decay in dileptons during the expansion. 

5 Medium modifications 

One of the most interesting aspects of studying low mass dilepton pairs is the possibility 
to learn something about in medium modifications of hadrons. One possibility is to see 
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the possible dropping of the mass of vector mesons |2], |7j as proposed by Pisarski [O 



as well as Brown and Rho [16]. In the context of the CERES data, this possibility has 



been studied in references [TT| as well as [I2~|. However, as far as the restoration of chiral 
symmetry is concerned, the dropping of the vector meson masses are at best an indirect 
signature. To leading order in the temperature chiral symmetry predicts no change in 



the mass of the vector mesons ]37|]. This finding is also confirmed by calculations using 
effective chiral Lagrangians [^, |39fl . 

In this article we want to concentrate on two medium effect which have been dis- 
cussed in detail in ref f|, |, ||, namely the reduction of the pion form factor due to the 
onset of chiral restoration and the effect of the in medium pion dispersion relation. As 
demonstrated in ref. || both effects combine to flatten the invariant mass distribution 
of dileptons in qualitative agreement with the CERES data. In ref. |5| the contributions 
of other production channels have not been taken into account. Also the effect of the 
experimental cuts have been ignored thus not allowing for a quantitative comparison with 
the data. Here we will attempt a first quantitative comparison with the data by including 
the above medium effects into our transport model. This is done in the following way: 
For each collision we determine the local temperature from the local covariant pion 
density 

where f{x,p) is the pion phase space distribution. Assuming local thermal equilibrium 
we obtain the temperature from 



/p(r)(2 



7T 



|2 



T(r) = f K 7 g ' (21) 

where p(r) denotes the pion density at position r. Given the local temperature we multiply 
the pion annihilation cross section by the ratio between in medium and free production 
as determined in ref. || || (see fig [l]). 

This procedure involves several approximations: First, we ignore any non equilibrium 
effects on the in medium corrections. This, however, may not be such a bad approxima- 
tion, since to leading order in the temperature the in medium corrections are really effects 
of the pion density, which we use to extract the temperature. Secondly, we do not prop- 
agate the pions according to their in medium dispersion relation, but only take its effect 
into account when determining the annihilation cross section as described above. This, 
however, leads only to a slight difference in the overall reaction dynamics as demonstrated 
in ref. |5(J . Finally, our largest and so far least controlled approximation is to assume that 



the pion-form factor does not depend on its momentum with respect to the heat bath. 
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Figure 1: Correction factor used in in medium calculation for the pion annihilation cross 
section 

The determination of the momentum dependence represents a rather involved computa- 
tion which we are currently trying to work out. However, considering the comparatively 
small contribution from the pion-pion annihilation (see below), we do not expect that 
the correction from the true momentum dependent form factor will alter our conclusions 
substantially. 

6 Results and Discussion 

Before we discuss our result in detail let us summarize our findings. By playing with 
the initial conditions for the hadronic system we could generate variations in the dilepton 
spectra, subject to the CERES acceptance cuts by at best a factor of two. A large portion 
of these variations arises from the contribution of the eta - Dalitz decay, and, therefore, 
is due to a principally measurable quantity. Considering the the large systematic and 
statistical errors of the present CERES data, we unfortunately cannot conclude that these 
present data put strong additional constraints on the hadronic equilibration dynamics. 
With some initial conditions (not including medium modifications) we are able to come 
within the lower end of the sum of statistical and systematic error of the CERES data and 
at the same time obtain a reasonable agreement with the HELIOS data. We furthermore 
find that the in medium modification considered give only an insignificant enhancement 
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of the total dilepton yield. The reason for that is that the pion annihilation, although 
important, contributes at best a third of the total yield in the relevant low mass region. 
Therefore a factor of two enhancement reduces to less than 30 % enhancement in the total 
dilepton spectrum. 



6.1 Dilptons from different inital states 

Let us now turn to the detailed discussion of our results. We have investigated many 
different initial conditions and the most representative of those are listed in table ([!]). 
Following the discussion in section |^ the characterization of the initial hadronic system 
can be divided into two major groups 

1. Momentum and configuration space (phase space). This includes choices for the 
formation time and possible initial radial flow. 

2. Relative abundances of particles 

The initial momentum space distribution is controled by the choice of rapidty distri- 
bution, transverse momentum distribution, initial flow velocity and, at least for the pions, 
a possible chemical potential. The only parameter we used to change the initial config- 
uration space is the longitudinal extend of the fireball, assume that initially transverse 
motion is small. Finally, we can generate possible configuration space and momentum 
space correlations by either chosing a finite formation time or, explicitly, by a finite value 
of Yi in eq. (0). 

We find that an initial number of about 600 pions (possibly 'hidden' in p etc mesons) 
is needed. The relative abundances are either taken as those found in nucleon-nucleon 
experiments (indicated by 'pp' in table ([[])) or as given by chemical equilibrium (indicated 
by ' chem 1 in table (|l|)). Whenever we have chosen chemical equilibrium abundances we 
also have assumed a vanishing formation time together with z-y correlations (Yi ^ 0) and 
a slightly more extended fireball (typically Zi = 2 fm, except sets 5 and 6) 

In figures ([3] - |8|) we show the resulting dilepton invariant mass spectra using the 
parameter sets given in table (p. To demonstrate the typical agreement we tried to 
achieve with the pion spectra and the proton and pion rapidity distribution the initial 
and final distributions for parameter sets 1 and 2 are shown in fig. (||]). For the pt proton 
spectra the slope parameter of the final spectrum (Tf ina i(prot) is given in table ([!]). These 



should be compared with an experimental slope parameter of about 2A0MeV [52]. Since 
the protons only indirectly affect the dilepton production we did not put too much an 
effort in fine tuning the resulting proton spectra. 

Parameter set 1 assumes that at the moment of hadronization the longitudinal position 
and the rapidity are correlated (YJ 7^ 0, see eq. (|T1|) ). This is similar to the initial 
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Set 1 


Set 2 


Set 3 


Set 4 


Set 5 


Set 6 


Set 7 


Set 8 


ratios 


P-P 


P-P 


P-P 


therm 


therm 


therm 


P-P 


therm 




350 


350 


350 


114 


156 


270 


1200 


540 


N P 


70 


70 


70 


96 


89 


95 


240 


307 




24 


24 


24 


60 


62 


31 


80 


214 


A T 

N ai 


12 


12 


12 


36 


25 


15 


40 


86 


r [rm/cj 





1 


0.5 











0.5 





T. ■ \Mc\f] 
J-imt [IVievj 




zuu 


zuu 


zuu 


i on 


1 fin 

1UU 


1 fin 

1UU 


1 on 
ryu 


[MeV] 











130 


130 








130 




1.6 








1.6 


1.6 


1.5 





1.6 


a 


0.55 


1.5 


1.5 


0.55 


0.6 


0.7 


1.8 


0.6 


Z l [fm] 


2 


1.5 


1.5 


2 


2.9 


22 


2.0 


3.4 


Yi(nucl) 


1.7 








1.7 


1.8 


1.8 





1.8 


a{nucl) 


0.6 


1.6 


1.6 


0.6 


0.7 


0.7 


1.9 


0.7 


Vflovi 














0.35 


0.4 





0.4 


T /w (prot.) [MeV] 


235 


230 


250 


240 


260 


250 


250 


290 



Table 1: Parameter sets used for dilepton calculations for S + Au. The line labeled 'ratios' 
indicates whether the relative abundances are taken according to proton proton data (p- 
p) or according to thermal equilibrium (therm). The last row gives the slope parameter 
of the final proton spectrum. 



conditions used by Li et al. [ll|. The relative particle abundances are assumed to be 
the same as in proton proton collisions. This leads to about initial 350 pions. With 
the longitudinal extent of Ri = 2Zi = 4fm and a radius of 3.5 fm the resulting initial 
pion density is then as large as pi n a = 2.3 fm -3 , which is certainly larger then one would 
consider reasonable for a transport descriptionP]. The initial density is reduced if one uses 
thermal particle ratios, which, for the same longitudinal extend, corresponds to parameter 
set 4. There we have assumed a fairly large chemical potential of /x w = 130 MeV. In this 
case the pion density reduces to = 0.75 fm -3 , which is an acceptable value for a 
transport description. Notice that the thermal pion density with that chemical potential 
and the temperature of T init = 260 MeV would be about twice as high. In this sense 
parameter set 4 also does not correspond to initial chemical and thermal equilibrium. 

Initial densities consistent with the thermal ones are considered in parameter sets 5 
and 6. Here we also have allowed for an initial transverse flow resulting in a lower initial 
temperature. The longitudinal extension of the fireball has then been fixed by requiring 



3 Notice, that due to the initial z — y correlations the effective pion density is, however, somewhat 
lower 
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Figure 2: Pion transverse momentum spectrum (left) and pion and proton rapidity dis- 
tribution (right) for parameter sets 1 and 2. The experimental data for the transverse 

(full circles) and the WA80 



momentum spectrum are from the NA35 collaboration 
collaboration |27j (open circles, errorbars left off). The data for the proton and pion rapid- 
ity distributions are from the NA35 collaboration The open squares are a reflection 
of the datapoints with Y > 2.6 with respect to Y = 2.6. 



that the initial hadron densities agree with thermal values for the respective chemical 
potentials. In these cases the initial pion densities are lower, due to the smaller initial 
temperature. They are pi n u = 0.68 fm -3 and pi n i t = 0.16 fm -3 for set 5 and 6, respectively 
A complementary, and maybe more realistic picture of the initial hadronic state, is to 
assume that the creation of a hadron from string fragmentation requires a certain time, the 
so called formation time (see section [|). The introduction of a formation time considerably 
reduces the initial effective hadron density because only hadrons which have been formed 
are allowed to interact. These are essentially only the particles in the same rapidity 
interval. Parameter sets 2 and 3 have a formation time of r = lfm/c and r = 0.5fm/c, 
respectively. In this case we did not give the particles an initial z — y correlation, since 
our physical picture is that particles are produced independently over the entire volume 
by the first hard nucleon-nucleon collisions. However, due to the finite formation time, for 
particles which are allowed to interact longitudinal position and rapidity are correlated. 
For the longitudinal extent we chose a size of Z\ = 1.5 fm, which corresponds to the 
Lorentz contracted longitudinal size of the combined sulfur and gold nuclei in the rest 
frame of the fireball. The initial temperature needed to reproduce the pion spectra is 
then T init = 200 MeV. 
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Finally we should point out that in all parameter sets except sets 5 and 6 the relative 
width of the rapidity distributions (see eq. flTTD ) between tt, p and u has been taken from 
proton proton data P5| . These ratios are Opjo^ = o^jo^ = 0.66. For the d\ we assume a 



slightly smaller ratio, namely cr^/ov = 0.5. 
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M [ GeV ] M [ GeV ] 



Figure 3: Dilepton production for parameter set 1 

Parameter sets 1 and 2 (see figures |3] and ^) give essentially the same results for the 
CERES dilepton spectrum. The differences seen in the comparison with the HELIOS 
data around the u mass demonstrates the effect of the formation time. In set 2, where 
we have a finite formation time of r = 1 fm the contribution of the pion annihilation is 
reduced for the forward rapidities where the HELIOS acceptance is located. 

The comparison with the calculations seems to suggest that the double hump structure 
in the HELIOS data is mainly due to rj Dalitz decays at low masses^ and the p / uj decays 
as well as pion annihilation at high masses. Consequently, possible small discrepancies at 

4 Notice, that HELIOS is measuring dimuons so that the minimum invariant mass is 200 MeV and, 
therefore, the rise in the rj Dalitz spectrum does not appear as in the CERES spectrum. 
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Figure 4: Dilepton production for parameter set 2 



low invariant masses (at the '^-hump') should only be taken seriously once a precise mea- 
surement of the rj abundance at forward rapidities is available. Without this information 
the disagreements seen in our results can always be 'fixed' by fine tuning the forward pion 
rapidity distribution and the hard end of the pion transverse momentum spectrum, since 
we use m t scaling to obtain the 77 mesons from the final pion distribution. 

Reducing the formation time to r = 0.5fm/c (Set 3, figure |5] ) increases the contri- 
bution from the pion annihilation, bringing our results close to the lower end of the sum 
of statistical and systematic error bars of the CERES data in the interesting mass region 
around 400 — 500 MeV. The agreement with the HELIOS data is pretty reasonable in this 
case except around an invariant mass of about 500 MeV, where we are about a factor of 
two too low, even if we fine tuned the 77 Dalitz contribution. In this case also the contri- 
bution of the 77-Dalitz is slightly larger then in Set 1. Because of the shorter formation 
time more flow is generated which in turn gives a slightly harder pion spectrum (which, 
however, is still in very good agreement with the data) and thus a few (20 %) more 77 
mesons. 
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Figure 5: Dilepton production for parameter set 3 

The one single process which contributes directly to the interesting region around 
M ~ 400 — 500 MeV is the Dalitz decay of the u meson. Since the creation cross section 
of uj meson via n-p fusion is rather small the only way to increase the number of u mesons 
is by assuming thermal particle ratios with a large pion chemical potential. This is done 
in parameter set 4, where we have assumed a pion chemical potential p w = 130 MeV. 
Comparing with parameter set 1, this clearly increases the dilepton yield in the relevant 
mass region and the main contribution is from the w-Dalitz decays. At the same time also 
the contribution of the direct decay of the u is increased, leading to an overprediction of 
the HELIOS data in the p-u mass region. 

As already mentioned, for parameter sets 5 and 6 we have adjusted the longitudinal 
size of the fireball such that the particle density agrees with the thermal one. We have 
also allowed for initial radial flow and, as a result, the initial temperature is smaller. Fur- 
thermore, unlike in the previous cases, the widths of the rapidity distributions have been 
taken to be the same for all mesons. Assuming thermal densities, although apparently 
more consistent, reduces the hadron dynamics to some extent, in particular in case of set 
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Figure 6: Dilepton production for parameter set 4 

6, where we have a longitudinal size of 44 fm. These are essentially freeze out conditions 
and it appears rather unlikely that hadrons do not interact and/or have expanded in the 
radial direction before that. In case of set 5 the longitudinal extent is not much larger 
than in the comparable set 1 and this may be a more realistic scenario, where hadrons are 
created in thermal equilibrium. Again, because of the large pion chemical potential in set 
5, there is a large fraction of initial omegas. Consequently, we find a sizeable contribution 
in the interesting mass region and an overshooting in the p uj region in case of the HELIOS 
data. 

Additional information may be obtained by studying transverse momentum spectra 
of the dileptons once data of sufficient statistics are available. In figure |9] we show the 
resulting transverse momentum spectra using the parameter set 3 for the invariant mass 
regions 400 MeV < M < 450 MeV (a) and 750 MeV < M < 800 MeV (b), respectively. We 
do not observe any significant structure in the spectrum and all contributions look more 
or less exponential. Only in the p - uj mass range (b) the contribution of the omega decay 
shows a significantly larger slope parameter ('temperature'). This is essentially a result 
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Figure 7: Dilepton production for parameter set 5 



of the transverse flow which is generated during the expansion. This affects mostly the u 
mesons, which decay predominantly after freeze out, i.e. after the transverse flow has been 
build up. Pion annihilation, the other strong contribution in this mass range, on the other 
hand, contributes mostly at the initial stage of the the expansion, and thus is not affected 
by flow. Therefore, already with the present mass resolution of the the CERES detector, 
the relative importance of the omega decays can be determined by a careful measurement 
of the slope parameter of the dilepton transverse momentum spectrum. This will then 
provide a crucial constraint to available model calculations. 

To conclude this part, it appears that without any in medium effects we are able to 
more or less reproduce, within errors, the CERES data as well as the HELIOS data. 
Except for the first three data points, our 'best' result, set 3, seems to be consistent with 
the CERES data, if the points are shifted towards the lower end of the systematic error. 
Looking at the HELIOS and CERES data together, there seems to be some enhancement 
at an invariant mass slightly above 500 MeV. However, given the large systematic errors 
of the CERES data, and the fact that the systematic errors of the HELIOS data are 
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Figure 8: Dilepton production for parameter set 6 



not given, it seems to be too early to declare the presence of any in medium effects. 
If, of course, an improved measurement would confirm the central values of the CERES 
data, from our investigations, we do not expect that any non-equilibrium effects of the 
hadronic phase would give that many dileptons in the relevant invariant mass region 
around 400 - 500 MeV. 

We have demonstrated that different initial hadronic configurations can lead to dif- 
ferences in the dilepton yield while still consistent with the hadronic observables. It also 
seems that the HELIOS data rule out initial partical ratios derived from chemical equi- 
librium based on a large pion chemical potential. Unlike the p, the decay of the u into 
dileptons is not expected to be reduced as a result of the resoration of chiral symmetry 
f44| . Also, a possible collisional broadening of the omega would not change this, because 



the experimental resolution is already much larger than any expected in medium width 
of the omega |]4"5| . 

From all the many configuration we have investigated (also those not reported here) it 
appears, however, that it is very unlikely to go above the lower end of the combined error 
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Figure 9: Transverse momentum spectra for dileptons using the CERES acceptance and 
parameter set 3. The ranges in invariant mass are 400 MeV < M < 450 MeV (a) and 
750 MeV < M < 800 MeV (b) 



of the CERES data the mass region around 400 — 500 MeV. As we shall discuss in the 
following section, we also do not expect that our proposed in medium correction would 
improve the situation considerably. 

6.2 Medium effects 

Let us now turn to the results including the proposed in medium modification of the 
pion annihilation cross section (see section ||). Since our in medium modification only 
affects the pion annihilation let us consider the parameter set 3, which has the biggest 
contribution from this channel. This then gives an upper limit of what we should expect 
for the different scenarios considered here. The resulting dilepton spectra are shown in 
figure together with the results without in medium effects, c.f. figure [5]. We only show 
the total yield and the contribution from the pion annihilation. All other contributions 
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are identical to those in figure |5|, since the way we implement the in medium corrections 
does not affect the expansion dynamics. While there is an enhancement in the pion 
annihilation contribution at low invariant masses and a suppression around the p-u mass, 
the effect in the total spectrum is hardly visible, especially in the interesting mass region. 
This is simply due to the fact that the pion annihilation contributes less than 1/3 to the 
total yield in this region and even an enhancement of a factor of two would increase the 
total spectrum by less than 30 %. 




M [ GeV ] M [ GeV ] 

Figure 10: Dilepton production for parameter set 3 with (full lines) and without (dashed 
lines) in medium modification of the pion annihilation cross section. 

Comparing with the changes in the dilepton yield we could generate by altering the 
initial conditions for the hadronic phase, the effect of the in medium corrections is small. 
We, therefore, have to conclude that it is probably impossible to learn something about 
the in medium pion dispersion relation from SPS-energy dilepton experiments unless the 
dynamics of the hadronic phase is extremely well known, and of course, much more precise 
measurements are feasible. 
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6.3 Prediction for Pb+Pb 



Currently, the system Pb + Pb at 150 GeV/A is investigated by the CERES collaboration 
so that it seems useful to try to make some predictions on what we would expect for 
the dilepton yield of this large system. Within our approach, this can only be a rough 
estimate, since final hadronic data for this system are not yet available. However, prelim- 
inary results for the rapidity distribution of negatively charged particles from the NA49 
collaboration |^6| as well as preliminary spectra from the NA44 collaboration are 
available. NA44 reports that the slope parameter for pion transverse momentum spec- 
trum changes only little when going from S+Pb to Pb+Pb. We, therefore, require that 
the final pion spectrum agrees with the S+Au data from the NA35 collaboration up to 
an overall constant. Our resulting rapidity distributions are shown in figure [11] together 



with preliminary data from the NA49 collaboration [46 




Figure 11: Pion and proton rapidity distribution for Pb+Pb. Data points are preliminary 
results for negatively charged hadrons from the NA49 collaboration. 

In figure [12|(a) we show our prediction for the Pb+Pb result, based on the CERES 
acceptance and normalization, i.e. dileptons per charged particles. The parameterizations 
we have used are given in table [I] (sets 7 and 8). They follow the same philosophy as sets 
3 and 5 for the S+Au case, respectively. In order to model a hadronic system generated 
in Pb+Pb collisions we increased the transverse radius to 6.5 fm and the number of total 
pions (~ 2000) and baryons (400). The relative abundances are the same as in sets 3 and 
5, respectively. 
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For a comparison, we also show the results for S+Au obtained with parameter set 3. 
We find, that if normalized to the number of charged particles, the dilepton spectrum 
hardly changes by going from S+Au to the heavier system of Pb+Pb. This result may 
be surprising at first sight, because in the heavier system many more pions are created 
which should give rise to more dileptons from pion annihilation. However, one has to be a 
little careful with this argument. The rate of dilepton production from pion annihilation 
is roughly proportional to the square of the pion density 

~Vt~ ~ Pn (22) 

where V is the the volume and r is the lifetime of the fireball and p n the average pion 
density. Consequently, the total number of dileptons produced divided by the number of 
charged particles (mostly pions) is linearly proportional the average pion density 

TP* {26) 



^charged ^charged 

Here we assumed that N charged ~ N^. It is probably fair to assume that the volume of 
the fireball scales like the volume of the initial nuclei and, therefore, more or less like 
A, the number of nucleons. But also the number of produced pions (negatively charged 
hadrons) seems to scale the same way. NA35 reports about 100 negative pion produced 
in S+S collisions f|7| and NA49 a preliminary number of about 700 in Pb+Pb collisions 



46| . Consequently, since both the volume and the number of pions seem to scale more 
or less with the size of the system, the average pion density reached in these collisions is 
essentially independent of the system size. Remains the lifetime of the fireball, r, which is 
probably larger for a bigger system but depends on the detailed expansion dynamics. The 
effect of the lifetime is demonstrated in figure |i~2|(b), where the contributions of the pion 
annihilation is plotted. If we compare the results for S+Au (set 5) with those for Pb+Pb 
(set 8), we find about a factor of 2 increase around the p-uj mass. In both calculations the 
initial pion density is roughly the same, so that difference comes from the longer lifetime 
of the larger system. 

Notice, that the CERES experiments accepts dileptons in the rapidity interval 2.1 < 

Y < 2.65 and, thus, favors the S+Au system, which has a fireball rapidity of about 

Y = 2.6, whereas Pb+Pb is located at Y = 3. This difference in acceptance is also the 
reason why the increase in the pion annihilation contribution is not reflected in the total 
yield, especially around the p-u mass. Dileptons from u and initial p mesons, which scale 
like the number of charged particles, contribute less. 

Again, our proposed in medium modification due to the pion dispersion relation leads 
only to a small correction (see dotted curve in figure |T2|(a) ). 
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Figure 12: Prediction for Pb+Pb. (a) Total yield, (b) contribution from pion annihilation 



To summarize our prediction, if there are no in medium effects such as dropping p 
masses, we expect that the dilepton spectrum for Pb+Pb subject to the CERES ac- 
ceptance and normalization should not differ significantly from that for S+Au. If the 
statistical and systematic errors are smaller for the Pb+Pb data, as expected 11481 , we 
would actually expect that in the mass region of 400 — 500 MeV the central values should 
come down somewhat as compared to the published S+Au data. If, on the other hand, 
the Pb+Pb data show a significant increase over S+Au, some additional, prehadronic 
production mechanism or strong in medium modifications such as dropping vector meson 
masses etc. would be called for. 



6.4 Baryon decays 

In our calculation we have not taken into account the possible Dalitz decays of any 
baryons. While we have taken into account the Delta resonance for the expansion dy- 
namics we have not included its Dalitz decay mainly because it will only contribute to 
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invariant masses lower than 300 MeV. However, at the temperatures considered also heav- 
ier baryon resonances will be present. For a temperature of 180 MeV we find that about 12 
% of the baryon number will be in states heavier than the Roper resonance (A r *(1440)), 
such as the JV*(1520), JV*(1535), the A(1600), and the A(1620). While the iV*(1440) 
hardly couples to photons (the branching ratio is smaller than a tenth of a percent) the 
higher nucleon and Delta states have considerable photon decay width. Also their masses 
are large enough such that the Dalitz decays could contribute to the relevant region around 
400 MeV. Since the lifetime or these resonances is rather short (r > 150 MeV) the domi- 
nant contribution to the dilepton spectrum will come during the lifetime r of the fireball, 
where we may assume a constant number of these resonances. The total contribution to 
the dilepton yield, therefore, is 

y\baryons 

Ajbaryons _ at j^baryons , *j 1 e+e~ (c)a\ 

iv e+e- ' 1 ^baryons L e + e - T ^'"baryons ^baryons \ > 

J- tot 

The second term takes into account the contribution after freeze. It can be ignored for 
our estimate since the lifetime of the resonances l/T^ yons is much smaller that that of 
the fireball (t). 

To give an upper limit for the contribution we assume a lifetime of r = 10 fm/c for the 
fireball and that the heavy baryons constitute 25 % of the total baryon number (which 
is a factor of two more than chemical equilibrium would suggest at a temperature of 
180 MeV). For the Dalitz decay width T b ^° ns we use the formula (4.8) of Wolf et al. |4j, 



which has been derived for the Dalitz decay of the Delta of a given mass M. The coupling 
we adjust by assuming that the partial width into photon plus nucleon is 1 MeV, which 
again is on the large side. For the mass we chose M = 1.6 GeV which is in between the 
nucleon and Delta excited states under consideration. Finally, as done in all calculations, 
we assume that there are 60 baryons in the fireball, so that we are dealing with 15 high 
mass baryons. 

As a reference we use the Dalitz decay of the u meson, since its contribution can also 
be simply estimated. Due to its small decay width, most u mesons decay after freeze out 
and, therefore, their contribution is given by 

iV e V = Njfgz. (25) 

1 tot 

In fig. ( [IB] ) we show the ratio between Dalitz decays from the baryon and that of the 
oj mesons, assuming that there are 24 uj mesons at freeze out. This number is comparable 
with our calculation presented in fig (|3]). We see that even with our optimistic assumption 
the contribution of the baryons to the dilepton spectrum is less the half of that from the 
Dalitz decay of the ui. It is probably fair to assume that the experimental acceptance cuts 
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apply equally to both. By comparing with the result shown in fig. (^) we thus conclude 
that the contribution of the baryon Dalitz decays only lead to a small correction ~ 10% 
to the results presented above, which may be slightly larger but certainly not significant 
in case of Pb+Pb. 
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Figure 13: Ratio of baryon-Dalitz decay over cu-Dalitz decay 



7 Conclusion 

In this article we have studied dilepton production form a hadronic system believed to 
be created in an SPS energy heavy ion collisions. We have addressed the question to 
which extent dilepton data provide information on the equilibration dynamics in the 
hadronic phase. To this end we have investigated several different initial conditions for 
the hadronic phase, subject to the constraint that the final hadronic observable agree with 
the data. The comparison with the dilepton data could rule out an initial state where 
the relative abundances of hadrons is given by chemical equilibrium based on a large 
pion chemical potential. Variations of the initial momentum and configuration space 
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distributions affected the dilepton invariant mass spectrum by not more than a factor of 
two, which is small compared to the error bars of the data. Therefore, with the presently 
available data, the question of the equilibration dynamics can not be settled. 

We were not able to find a scenario which would give enough dilepton in the mass 
region around 400 — 500 MeV to reach the central values of the CERES data. Our 'best' 
results agrees with the lower end of the sum of statistical and systematic errors in the 
region. Considering the large experimental uncertainties we can not conclude that any in 
medium modifications are needed in order to explain the present data. 

We predict essentially the same dilepton spectrum for the larger system Pb+Pb as for 
S+Au, provided the CERES acceptance and normalization are used. Since for this system 
smaller experimental errors are expected, these data may give very valuable information 
about possible in medium corrections such as dropping vector meson masses. 

We have investigated the effect of an in medium pion dispersion relation and pion 
form factor on the dilepton spectrum. While the effect is visible in the pion annihilation 
channel, it hardly affects the total dilepton yield. Therefore, unfortunately it seems rather 
unlikely that the presence of an in medium pion dispersion relation can be detected in 
SPS-energy heavy ion collisions. 

Our conclusions could be much more sharpened if precise measurements of the the 77 
mesons and the u mesons would be available. The latter, of course, could be obtained best 
with a considerably improved mass resolution of the dilepton spectrometer. However, a 
careful measurement of the transverse momentum spectrum in the p - 00 mass region would 
already provide some information on the relative importance of the omega decays. Further 
constraints could be obtained by looking at pion correlation (HBT) data. These may 
provide valuable information about the space-momentum correlations and thus restrict 
our choices of initial conditions. However, considering, the rather weak sensitivity of the 
dilepton data with respect to the variation of the initial phase space, we do not expect 
much additional insight from such an investigation. 
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